Dependence of Striped Domain Growth on Background Wavenumber 
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We report on the growth of domains of standing waves in electroconvection in a nematic liquid 
crystal. An ac voltage is applied to the system, forming an initial state that consists of travelling 
striped patterns with two different wavenumbers, zig and zag rolls. The standing waves are gen- 
erated by suddenly applying a periodic modulation of the amplitude of the applied voltage that is 
approximately resonant with the travelling frequency of the pattern. By varying the modulation 
frequency, we are able to vary the steady-state, average wavenumber. We characterize the growth 
by measuring the structure factor for the pattern, the amplitude of the pattern, and the domain- 
wall length. An interesting feature of the growth is that the evolution of the domain-wall length 
is strongly affected by changing the selected wavenumber; whereas, the evolution of the structure 
factor is less sensitive to such changes. 

PACS numbers: 47.54.+r,64.60.Cn 
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I. INTRODUCTION 

When a system experiences a sufficiently large, sud- 
den change of an external parameter, referred to as a 
quench, the current state of the system will lose stabil- 
ity and a transition will occur to a new state. When 
there is more than one possibility for the new state (such 
as spin up or down in a magnetic system), domains will 
form. At late times, these domains coarsen, or grow. 
The dynamics of this coarsening process has been the 
subject of substantial research and is important in many 
applications Q. In general, the domains are character- 
ized by a single length scale, l(t), that grows according 
to a power-law, lit) ~ t n . Here n is referred to as the 
growth exponent. This process is usually referred to as 
domain coarsening. Some examples of systems for which 
domain coarsening has been studied in detail 
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theoretically or experimentally, include the Ising model, 
binary fluids, liquid crystals, diblock copolymers 0,0. 
and electroconvection Q. 

There are a number of ways to categorize systems un- 
dergoing coarsening. First, one must distinguish between 
cases where the order parameter is conserved, such as in a 
binary fluid, and where it is not conserved, such as in the 
Ising model. Second, the domains that are growing may 
be uniform, such as in both the Ising model and binary 
fluids, or they may consist of a pattern, such as stripes in 
the case of diblock copolymers and electroconvection. Fi- 
nally, the initial and final states can either be thermody- 
namic equilibrium states or steady states of the system. 
Most of the studies have focused on the case where the 
initial state prior to the quench and the final state well 
after the quench are thermodynamic equilibrium states 
[]}■ We will refer to this type of system as a thermo- 
dynamic system. For isotropic, thermodynamic systems, 
where the domains are uniform, the exponents are well 
known. For systems with a non-conserved order parame- 
ter, n — 1/2. Systems with a conserved order parameter 
have n = 1/3. In this report, we focus on coarsening in 
a system in which the domains are composed of stripes. 



An added feature of the system we consider is that the 
initial and final states are not states of thermodynamic 
equilibrium, though they are steady-states of the system. 
We refer to this as a driven system. 

For systems in which the domains are composed of 
stripes, theoretical studies of the power-law scaling of 
the late-time domain growth has focused on simulations 
and analysis of the Swift-Hohenberg equation 0, IE S IS 
0, fToL ITl| . Experiments have been performed in both a 
thermodynamic system, diblock copolymers 2. 3], and a 
driven system, electroconvection 0, ^3 ■ This work has 
raised two related issues. First, two length scales are 
important in the problem: the wavelength of the pat- 
tern and the scale of a domain. This has an interesting 
impact on the measured growth exponents. Traditional 
studies of domain growth focus on the structure factor, 
S(k). This measures both variations in wavenumber and 
domain size, with information about both lengths scales 
mixing in a nontrivial way. For striped systems, one can 
also consider orientational correlation functions, defect 
densities, or domain wall lengths. These quantities are 
independent of the local wavenumber and are expected 
to provide information about domain size only. 

Numerical simulations have considered both the scal- 
ing of the structure factor, S(k), of the system and the 
orientational correlation function. The evolution of the 
structure factor is found to be consistent with a length 
scale that grows as i 1 / 5 for both potential and nonpo- 
tential versions of the Swift-Hohenberg equation [(J |7j • 
However, the growth exponent determined from the ori- 
entational correlation function is 1/4 for the potential 
case and 1/2 for the nonpotential case |(|. The orienta- 
tional correlation function and defect densities are mea- 
sured experimentally in diblock copolymers Q. The re- 
sults agree with the potential case of the Swift-Hohenberg 
model, n — 1/4. For a particular case of electroconvec- 
tion, it was shown that the wavenumber and domains 
evolve differently Q . It was found that the growth expo- 
nent measured using S(k) is consistent with n = 1/5 and 
measurement of the domain size gave n = 1/4. This is 
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in agreement with simulations for the potential case [6(. 
In addition, the distribution of local wavenumbers was 
measured separately. This distribution was essentially 
independent of time, providing an explanation for the 
difference between S(k) and domain wall measurements. 

The second issue associated with striped domains is the 
interaction between defects in the pattern and the back- 
ground stripe pattern. This has been studied using the 
Swift-Hohenberg equation. One important consequence 
is that the growth of domains can be strongly affected 
by pinning of grain boundaries 0, 0| . in these simula- 
tions, the measured growth exponents depends on the 
depth of the quench, something that does not occur in 
uniform systems. The growth exponents are also found 
to depend on the level of noise in the system 0, [^} . This 
is attributed to the fact that noise can provide a source 
of "depinning" for the defects. Finally, the evolution of 
the average wavenumber provides a possible explanation 
for differences in simulations of the potential and nonpo- 
tential versions of the Swift-Hohenberg equations. The 
average wavenumber in each case evolves to different val- 
ues, and this was suggested as one possible explanation 
for the two different growth exponents . 

We report on an experimental study of the impact of 
changes in the average wavenumber on the growth of do- 
mains. We utilize a particular example of electrocon- 
vection that is ideally suited for controlling the average 
wavenumber of the pattern: electroconvection in the ne- 
matic liquid crystal 152 [l%ll4| . A nematic liquid crystal 
is a fluid in which the molecules are aligned on aver- 
age [l^. The axis along which the average alignment 
points is referred to as the director. For electroconvec- 
tion [lfj, a nematic liquid crystal is placed between two 
glass plates that have been treated so that the director is 
everywhere parallel to the plates. An ac voltage is applied 
perpendicular to the plates. Above a critical value of the 
voltage, a pattern of stripes forms. Because the system 
is anisotropic (the director selects a preferred axis), the 
states are usually characterized by the angle between the 
wavevector of the pattern and the undistorted director, 
6. States with the same wavenumber and angle ±9 are 
degenerate and referred to as zig and zag, respectively. 
For the system discussed in this paper, the initial pat- 
tern consists of four travelling stripe states ^?|: right - 
and left-travelling zig and zag rolls. It has been estab- 
lished for this system that above a critical value for a 
modulation of the applied voltage at twice the travelling 
frequency, the system consists of either standing zig or 
zag rolls |14j. Therefore, a sudden quench of the system 
from below this critical value to above it results in the 
nucleation of domains of standing zig and zag rolls, which 
proceed to coarsen. 

Previous measurements of growth exponents for this 
system focused on quenches that used exactly twice the 
travelling frequency. Reference £| reports a value of 
n = 1/4 for the growth of the domain size and n — 1/5 for 
the scaling of S(k). A single initial experiment for a mod- 
ulation frequency slightly different from twice the trav- 



elling frequency suggested similar behavior. However, a 
systematic exploration of changes in the wavenumber, by 
vary ing the modulation frequency, is possible for this sys- 
tem |l4| . and has not been carried out. In this paper, we 
report on such a study. The rest of the paper is orga- 
nized as follows. Section II provides the details of the 
experimental system and the techniques used to analyze 
the domain growth. Section III presents the results of 
the experiment. Section IV provides a discussion of the 
results. 



II. EXPERIMENTAL DETAILS 

We used the liquid crystal 152 0] doped with iodine. 
Due to the nature of 152, 8% iodine by weight was used 
initially. However, a significant fraction of the iodine 
evaporates before the cells are filled, so the final concen- 
tration is not well controlled. We used commercial liquid 
crystal cells [lj| composed of two transparent glass plates 
coated with a transparent conductor, indium-tin oxide. 
The cells were 2.54 cm x 2.54 cm, with the conductor 
forming a square 1 cm x 1 cm in the center of the cell. 
The glass was also treated with a rubbed polymer to ob- 
tain uniform planar alignment of the director. We define 
the direction parallel to the undistorted director to be the 
x axis and the direction perpendicular to the undistorted 
director to be the y axis. 

The apparatus is described in detail in Ref. [2(|. It 
consists of an aluminum block with glass windows. The 
temperature of the block was kept constant to ±0.001 °C. 
Polarized light is applied from the bottom of the cell. 
Above the cell, there is a ccd camera to capture the image 
of patterns using the standard shadowgraph technique 

m 

We apply an ac voltage of the form V(t) = y/2[V + 
V m cos{uj rn t)]cos(uidt) perpendicular to the glass plates. 
Here V a is the amplitude of the applied voltage in the 
absence of modulation, V m is the modulation amplitude, 
u) m is the modulation frequency, and LOd is the driving 
frequency. For all of the experiments reported on here, 
U!d/2ir — 25 Hz. The critical voltage, V c is defined as the 
value of V for which the system makes a transition from 
a uniform state to the superposition of travelling rolls in 
the absence of modulation. There are two relevant control 
parameters: e = (V /V c ) 2 — 1 and b = V m /V c . 

For the experiments discussed here, the system is ini- 
tially at e = 0.04 and 6 = 0, i.e. above onset with no 
modulation. As mentioned in Sec. I, there are only four 
modes needed to describe the observed convection rolls: 
zig and zag, right and left travelling oblique rolls. This 
state is an example of spatiotemporal chaos in that the 
amplitude of the four modes varies irregularly in space 
and time |22|. For each e and u> m , there exists a crit- 
ical value of 6, b Cl above which the system consists of 
standing zig and zag rolls. A quench corresponds to a 
sudden change of b to some value b > b c . At and below 
V c , the modulation frequency with the strongest coupling 
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to the pattern, and correspondingly smallest value of b c , 
is twice the Hopf frequency, oj^. The Hopf frequency is 
the travelling wave frequency at V c . As one increases e, 
the natural frequency of the pattern shifts. Because we 
worked at a finite value of e, the smallest value of b c oc- 
curs for a modulation frequency, uj m — 2lu* , where to* is 
the natural frequency of the pattern at the given value 
of e. Because it is the relative shift in frequency that 
matters, we will consider the behavior as a function of 
Sf = f* - /m/2, where /* = uj*/2it and f m = u m /2n. 
Figure is a plot of b c as a function of Sf for e = 0.04. 
We focus almost exclusively on the negative values of Sf 
because for Sf > +0.01 Hz, the pattern stabilized by 
modulation is standing waves that are a superposition of 
zig and zag rolls, and not individual domains of standing 
zig and zag rolls. 

Quenches were performed for a range of Sf. The final 
value of b in each case was selected so that the quench 
depths all corresponded to Ab = 0.02, as measured from 
the onset of standing waves (solid curve in Fig. ^| . The 
camera was triggered to take the images at the same 
point in the modulation cycle. This ensured that the pat- 
tern was imaged at the same point in the standing wave 
cycle. Each image covered a region of 4.1 mm x 3.1 mm 
of the cell. Time is scaled by the director relaxation time, 
Td = 71 cP /(tt 2 Kh) — 0.2 s. Here 71 is a rotational vis- 
cosity; Ku is the splay elastic constant; and d is the 
thickness of the cell. 

An example of a small region of the pattern after a 
quench is given in Fig. [21 Two different quenches are 
illustrated: Sf = 0.008 Hz and Sf = -0.23 Hz. The im- 
ages for Sf = 0.008 Hz illustrate the type of patterns on 
which we focused: standing zig and zag rolls separated by 
domain walls. The domain walls are relatively "thick" , 
and the superposition of zig and zag rolls in a domain 
wall is visible. The quench at Sf = —0.23 Hz illustrates 
qualitatively different behavior. In this case, there are 
regions of superposition that are "domains" in their own 
right, though most of the regions of superposition can 
be considered to be domain walls. As we will show, this 
qualitative difference is apparent in various characteri- 
zations of the growth. As the focus of the paper is the 
impact of varying the wavenumber on the domain growth 
in the zig/zag system, we will only report on this single 
quench in the regime where domains of the superposition 
of zig and zag rolls occurs. Further increasing \Sf\ results 
in patterns that are entirely the superposition of zig and 
zag rolls. Future work will consider in more detail this 
other class of patterns. 

A number of features of the domains were analyzed. 
First, the power spectrum, S(k), of each image was 
computed. This was used to determine the average 
wavenumber (< k >= J kS(k)dk/ J S(k)dk) for the 
zig and zag rolls. The widths of the fundamental 
peak in S(k) in both the x (a x ) and y (a y ) direc- 
tions were determined using the second moment of S(k). 
This gives a x = yj < k% > — < k x > 2 , with < k? x >= 
J k^S(k)dk/ J S(k)dk, and a similar definition for a y . 



(Note: all of the integrals were taken as sums over a fi- 
nite region around the fundamental peak, as described in 
Ref. [J.) a x and a y contain information about the typi- 
cal size of a domain and the distribution of wavenumbers, 
in the respective directions. 

We also analyzed the images by first converting them 
to a scaled image of zig and zag domains, referred to as 
the angle map. This was done using the procedure de- 
tailed in Ref. [4j . Briefly, the local angle of the wavenum- 
ber is computed using the algorithm reported in Ref. • 
This is used to produce an image with black and white 
domains representing regions of zig and zag, respectively. 
Regions of superposition of zig and zag rolls (the domains 
wall for most of the quenches) are represented by gray. 
Figure [3] illustrates a typical image (Fig. |3Ja,) and the re- 
sults after conversion to an angle map (Fig. |3Jd). Once 
these images are produced, two measures of the domain 
growth that are independent of variations in the local 
wavenumber are computed. First, one can determine the 
total domain wall length per viewing area, L. This quan- 
tity scales as the inverse of a characteristic domain size 
(the area of the viewing scales as a length squared) . In 
our case, L was defined as the total number of gray pix- 
els divided by the total area of the image. For all of the 
quenches except the ones at Sf = —0.23 Hz, the width of 
the domain walls was constant during the domain growth. 
Therefore, this is a consistent measure of the total do- 
main wall length. For the quench at Sf = —0.23 Hz, we 
still used this definition for consistency, even though it 
represents a transitional behavior. In this case, some of 
the regions of superposition of zig and zag rolls (gray pix- 
els) are more appropriately considered domains in their 
own right. Still, most of the gray pixels corresponded to 
domain walls between zig and zag domains. We also com- 
puted the power spectrum of the angle map image, S(q). 
In this case, the peak is centered at q = 0, as the angle 
map domains are uniform. Therefore, the width of the 
peak is related to the typical domain size only. Again, we 
separately considered the width in the x direction (w x ) 
and y direction {w v ) 1 where the widths are defined in the 
same manner as for S{k). 



III. RESULTS 

Figure 0] shows the time evolution of the average 
wavenumber (< k >) for the different quenches of inter- 
est. A number of features of the evolution are evident. 
First, as claimed, the steady-state value of the < k > is 
dependent on Sf (see Fig.0^,). However, for all quenches, 
the < k > is still evolving even at late times. After 
t/Td pa 1000, the evolution is consistent with a very slow 
logarithmic growth: < k >= a \og w (t/Td), with a « 0.01. 
This is illustrated in Fig. However, for all but the 
largest value of \Sf\, there is a weak dependence on when 
the logarithmic growth begins (with values ranging from 
2.7 < log 10 (t/rd) < 3.1) and on the coefficient (a) of 
the logarithmic growth (with values ranging from 0.006 
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to 0.02. In general, for larger \6f\, evolution consistent 
with logarithmic time dependence occurs at a later time. 
As referred to in Sec. II, for Sf = —0.23 Hz, the pattern 
is qualitatively different. In this case, the evolution of 
< k > is also very different, with logarithmic growth set- 
ting in significantly later (log 10 (i/r<z) ~ 3.5) and with a 
coefficient of 0.05. 

Figure [5] shows the evolution of the width of the power 
spectrum for both the pattern (a x and o~ y ) and the an- 
gle map (w x and w y ). There are a number of features 
that become apparent in these plots. First, if one were 
to assume that a scaling regime is reached, the growth 
exponent measured from S(k) is very different from that 
measured for S(q). Based on a x and a y , one observes 
very slow growth for the system. For late times, assum- 
ing a ~ t n , the exponents n are on the order of 0.1. In 
contrast, the length scales characterized by the inverse of 
w x and w y grow faster. In this case, one obtains growth 
exponents of n = 0.33 or larger. This difference is consis- 
tent with earlier results comparing S(k), which includes 
information about the variation in local wavenumber, to 
more direct measures of the domain size |6j . 

The second feature of the growth is the differences be- 
tween the growth in the x- and y-directions. For the 
quench with Sf = 0.008 Hz, both w x and w y have the 
same magnitude and time dependence. This suggests 
that the system is evolving in a relatively isotropic man- 
ner. For o~ x and o~ y , the time evolution is similar, though 
the magnitudes are different. This suggests that there 
is more anisotropy in the local wavenumber distributions 
than there is in the domains. This anisotropy contributes 
to the difference in magnitude between a x and o~ v , but 
does not impact the time evolution. For the case of 
Sf = —0.23 Hz, the situation is roughly reversed. In 
this case, a x and o~ y have roughly the same magnitude 
throughout the evolution. In contrast, w x and w y have 
different magnitudes and slightly different time evolution, 
at late times. This qualitative difference is probably due 
to the increased area of the pattern that is given by a 
superposition of zig and zag rolls, and thus, the change 
in the nature of the domain boundaries. 

As an alternative measure of the domain growth, we 
also considered the total domain wall length present in a 
fixed viewing area. This measure is an isotropic measure 
of the domain growth. Its time evolution is a good mea- 
sure of how rapidly the overall size of a domain changes. 
However, it does not provide any information about the 
relative rate of change in the x- and y-direction. Figure|B| 
shows the evolution of the total domain wall length for a 
number of values of Sf. Also shown are two straight lines 
that provide a guide to the eye. The dashed line has a 
slope of —0.7, and the solid line has a slope of —0.3. It 
is clear from this plot that the evolution of the domain 
size occurs more rapidly for the larger values of \Sf\. 

At this point, it is important to emphasize one of the 
challenges of interpreting data used to characterize do- 
main growth. In studies of domain coarsening, it is 
standard to assume that a late-time scaling regime ex- 



ists in which the lengths characterizing the domain size 
grow as a power-law in time. One experimental diffi- 
culty with this assumption is determining when such a 
scaling regime is reached. This is especially true in the 
studies reported on here, where the average wavenumber 
is still evolving in time for all of the conditions that we 
studied (even if it is only logarithmic in time). However, 
independent of whether or not a true scaling regime has 
been reached, for similar times during the evolution of 
the domains, the data is consistent with a power-law of 
the form l(t) ~ t n , where I is some characteristic length. 
Therefore, the values of n provide information about the 
relative rate of growth of domains under different cir- 
cumstances. This is our main motivation for looking at 
n. Any claims of a true scaling regime would be pre- 
mature at this point, though the consistency of the data 
with power-law behavior suggests that it is possible such 
a regime has been reached. 

Figure summarizes the impact of variations in the 
wavenumber on the domain growth by comparing growth 
exponents n for w x , w y , and the total domain wall length 
per viewing area L. The results from a x and o~ y exhib- 
ited no strong dependence on Sf. One observes that the 
growth of the domains is significantly faster for larger 
| Sf |, as measured by the total domain wall length in a 
given area. In this case, n increases in magnitude mono- 
tonically with \Sf\. However, the trends in w x and w y 
are less obvious. If anything, one observes a small change 
in the magnitude of n for w x and uu y as \Sf\ is increased. 
Also, the growth appears to become more anisotropic as 
\Sf\ increases. 



IV. SUMMARY 

The experiments reported on here demonstrate the 
close connection between the dynamics of domain growth 
in patterned systems and the evolution of the background 
wavenumber. We performed a series of quenches in which 
the steady-state average wavenumber after the quench is 
varied. The resulting domain growth depends on the 
evolution and steady-state value of < k >. However, 
different measures of the domain growth are impacted to 
different degrees. A direct measure of the typical domain 
size, the average domain wall length in a fixed viewing 
area, showed the strongest dependence on the average 
wavenumber, with the rate of growth increasing as the 
wavenumber is varied away from its "natural" value. 

If one assumes that the average domain wall length 
scales as a power in time, exponents in the range of 0.3 to 
0.7 are observed. As discussed in the introduction, such 
a large variation in growth exponent is consistent with 
simulations of the two different Swift-Hohcnberg equa- 
tions 0. In Ref. 0, the calculated exponent is found to 
increase from 0.25 to 0.5 with a variation in wavenum- 
ber of approximately 20%. This difference in average 
wavenumber was suggested as the source of the differ- 
ent exponents. The work in Ref. |6| does not directly 
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apply to this system because it considered the isotropic 
Swift-Hohcnbcrg equation, and our system is anisotropic. 
However, the similarity between the two results suggests 
that it is worth further work to determine how universal 
the impact of wavenumber variation is on domain growth 
in patterned systems. 

An interesting feature of the growth is that as one 
forced the average wavenumber to be different from the 
"natural" wavenumber of the system, the rate of growth 
was increased. In our case, the maximum variation in 
wavenumber was only 6%, but the exponent changed 
from 0.3 to 0.7. If these results carry over to other pat- 
terned systems, such as diblock copolymers, one can have 
a relatively significant impact on the rate of growth of 



domains by inducing a change in the average wavenum- 
ber away from the "natural" wavenumber. This increase 
in the rate of domain coarsening could have a practical 
impact in various processing situations. 
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FIG. 1: Onset curve for standing waves at a value of e — 0.04 
as a function of Sf = /* — /m/2. Here /* is the natural 
travel frequency of the pattern at e = 0.04 and f m is the 
modulation frequency. The parameter b = V m /V c , where V m 
is the modulation voltage and V c is the critical voltage for 
the onset of a pattern in the absence of any modulation. The 
solid line is a fit to the data used to determine /*. 





FIG. 3: (a) Image taken 384 s after a quench with Sf = 
0.008 Hz. (b) The same image converted to an "angle" map. 
The black and white regions are zig and zag rolls, respectively. 
The gray regions are the boundaries between domains. The 
scale bar in (a) represents 0.4 mm. 



d ■ 













FIG. 2: Four images of domain growth. Images (a) and (b) 
are with a modulation frequency Sf = 0.008 Hz, at t — 96 s 
and t — 600 s after the quench, respectively. Images (c) and 
(d) are with a modulation frequency Sf = —0.23 Hz, at t = 
107 s and t = 667 s after the quench, respectively. The solid 
bar in (a) is 0.2 mm. 
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FIG. 4: (a) Average wavenumber for the zig standing rolls 
as a function of time after the quench. Time is scaled by 
the director relaxation time, ra- The symbols are as follows: 
(□) Sf = 0.008 Hz, (o) Sf = -0.012 Hz, (a) Sf = -0.022 Hz, 
(V) Sf = -0.072 Hz, and (0) Sf = -0.23 Hz. The wavenum- 
ber for the zag rolls shows similar behavior, (b) An expanded 
view of the average wavenumber for Sf = 0.008 Hz. The solid 
line is a guide to the eye illustrating the regime where the 
evolution of < k > is consistent with logarithmic growth in 
time. 
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FIG. 5: Width of the power spectrum as a function of time 
for two different values of Sf. Time is scaled by the director 
relaxation time r d . (a) Sf = 0.008 Hz (b) Sf = -0.23 Hz. 
In both cases, the open symbols are for the power spectrum 
computed on the image with the stripes, with o x given by the 
squares and a y given by the circles. The closed symbols are 
for the power spectrum of the image of the angle map, with 
w x given by the squares and w y given by the circles. 
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FIG. 6: Total domain wall length as a function of time 
for all the quenches studied. Time is scaled by the di- 
rector relaxation time, Td- The symbols are as follows: 
(□) Sf = 0.008 Hz, (o) Sf = -0.012 Hz, (a) Sf = -0.022 Hz, 
(V) Sf = -0.072 Hz, and (0) Sf = -0.23 Hz. The lines are 
provided as a guide to the eye, with the dashed line having a 
slope of —0.7 and the solid line having a slope of —0.3. 
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FIG. 7: Summary of measured scaling exponents (n) for both 
the domain wall length (solid squares), w x (open triangles), 
and w v (open circles). 



